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ABSTRACT 

Using a 3D MHD simulation, we model the quasi-static evolution and the 
onset of eruption of a coronal flux rope. The simulation begins with a twisted 
flux rope emerging at the lower boundary and pushing into a pre-existing coro- 
nal potential arcade fleld. At a chosen time the emergence is stopped with the 
lower boundary taken to be rigid. Then the coronal flux rope settles into a 
quasi-static rise phase during which an underlying, central sigmoid-shaped cur- 
rent layer forms along the so called hyperbolic flux tube (HFT), a generalization 
of the X-line conflguration. Reconnections in the dissipating current layer effec- 
tively add twisted flux to the flux rope and thus allow it to rise quasi-statically, 
even though the magnetic energy is decreasing as the system relaxes. We examine 
the thermal features produced by the current layer formation and the associated 
"tether-cutting" reconnections as a result of heating and field aligned thermal 
conduction. It is found that a central hot, low-density channel containing recon- 
nected, twisted flux threading under the flux rope axis forms on top of the central 
current layer. When viewed in the line of sight roughly aligned with the hot chan- 
nel (which is roughly along the neutral line), the central current layer appears 
as a high-density vertical column with upward extensions as a "U" shaped dense 
shell enclosing a central hot, low-density void. Such thermal features have been 
observed within coronal prominence cavities. Our MHD simulation suggests that 
they are the signatures of the development of the HFT topology and the associ- 
ated tether-cutting reconnections, and that the central void grows and rises with 
the reconnections, until the flux rope reaches the critical height for the onset of 
the torus instability and dynamic eruption ensues. 
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Introduction 



CMEs and eruptive flares are believed to result from a sudden, explosive release of 
the free magnetic energy stored in the previously quasi-equilibrium, twisted/sheared coronal 
magnetic field (see e.g. reviews by [Forbes et"aL][2006l pienl[20TT1 ) . Many CME models have 
considered a magnetic fiux rope containing helical field lines twisting about each other in the 



corona as the basic underlying magnetic field structure for CME precursors (e.g. Titov and 



Demoulin 1999; Low 2001; Gibson and Fan 2006). Recent observational studies of coronal 



prominence cavities support the picture of twisted coronal fiux ropes as CME precursors (e.g. 
Wang and Stenborg||2010 ). Theoretical studies have shown that a coronal fiux rope confined 
by an external potential magnetic field can become unstable to the torus instability and erupt 
when it reaches a critical height where the ambient potential field decreases with height too 



steeply (Kliem and Torok 2006; Isenberg and Forbes 2007; Demoulin and Aulanier 2010). 



MHD simulations have been carried out to study the critical condition and the non-linear 



evolution of the torus instability of various 3D coronal fiux rope structures (e.g. Torok and 



Kliem 2007; Fan and Gibson 2007; Aulanier et al. 2010; Fan 2010). 



Through a sequence of 3D MHD simulations of the evolution of the coronal magnetic 
field resulting from the emergence of a twisted magnetic fiux rope into a pre-existing coro- 



nal arcade field, Fan (2010) studied the conditions that lead to a dynamic eruption of the 



resulting coronal fiux rope. It is found that the critical condition for the onset of eruption 
is for the center of the fiux rope to reach a critical height at which the corresponding po- 
tential field declines with height at a sufiiciently steep rate, consistent with the onset of 
the torus instability of the fiux rope. The simulations show that after the fiux emergence 
is stopped, the coronal fiux rope first settles into a quasi-static rise phase with an under- 
lying, central sigmoid-shaped current layer developing. Reconnections in the current layer 



is found to effectively add twisted fiux to the fiux rope (as quantified by eq. (14) in Fan 



2010), allowing it to rise quasi-statically, even as the magnetic energy is declining. We thus 



call these reconnections "tether-cutting" reconnections in the sense that they are effectively 
reducing the anchoring of the fiux rope by dissipating the underlying, central current layer 
whose current attracts the volume current in the fiux rope (Low and Zhang 2002). When 



the fiux rope rises quasi-statically to the critical height for the onset of the torus instability, 
then dynamic eruption ensues. To identify the thermal signatures that may result from the 
formation of the current layer and the associated reconnections, we carry out an MHD sim- 



ulation with a similar setup as those in Fan (2010), but with a more sophisticated treatment 



of the thermodynamics. Instead of assuming a simple isothermal evolution, we solve the 
total energy equation in conservative form for a perfect gas with an adiabatic index 7=1.1. 
In this way the non-adiabatic heating due to the dissipation of kinetic and magnetic energies 
by the numerical diffusions are implicitly incorporated in the thermal energy. Also thermal 
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conduction along the magnetic field lines is included in the energy equation. We find from 
this new simulation that the main consequence of the tether-cutting reconnections in the 
current layers during the quasi-static phase is the formation and enlargement of a hot, low 
density channel, containing reconnected, twisted fiux threading under the axis of the fiux 
rope. When viewed above the limb along a line of sight that is roughly aligned with the 
channel, the thermal signature would manifest in SDO/AIA observations as a hot central 
void enclosed in a "U" shaped dense shell on top of a relatively dense structure inside the 
flux rope. This signature inside the flux rope is consistent with the recent observations of the 



thermal sub-structures inside coronal prominence cavities by Reeves et al. (2012); Berger 
( |2012| ); [Rifflier et al.| ( [20TT[ ). 

We organize the remainder of the paper as follows. In Section [2| we describe the MHD 
numerical model with the more sophisticated treatment of the thermodynamics compared 
to the previous simulations in Fan (2010). In Section [3] we describe the simulation result, 
focusing on the thermal signature that develops inside the flux rope and compare it with 
observations. We summarize and discuss the conclusions in Section [H 



2. Model Description 



For the simulation carried out in this study, we solve the following magneto- hydrodynamic 
equations in a spherical domain: 
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q = -/^oT^/^bb • VT (9) 

In the above v is the velocity field, B is the magnetic field, p, and T are respectively the 
plasma density, pressure and temperature, 5, e i?, /x, and 7, are respectively the internal 
energy density, the total energy density (internal+kinetic+magnetic), the gas constant, the 
mean molecular weight, and the adiabatic index of the perfect gas, q is the heat fiux due 
to field aligned thermal conduction, with the field direction denoted by the unit vector b 
and with /^o = 10~^ erg K"''/^ cm~^ s~^, and G and Mq denote the gravitational constant 
and the mass of the Sun. We have assumed a perfect gas with a low value of the adiabatic 
index: 7 = 1.1, which allows the coronal plasma to better maintain its (high) temperature 
during adiabatic expansion. No explicit viscosity and magnetic diflFusion are included in 
the momentum and the induction equations. However numerical dissipations are present at 
regions of sharp gradient, and since we are solving the total energy equation of the plasma in 
conservative form (eq. [5]), the non- adiabatic heating due numerical dissipation of kinetic, 
and magnetic energies is implicitly put back into the internal energy. We also explicitly 
include the non-adiabatic effect due to the field- aligned thermal conduction (eqs. [5] and 
[9]) in the energy equation, but in a limited way. Given the temperature and density in the 
solar corona, the field aligned thermal conduction would pose a severely limit on the time 
step of numerical integration. Here we artificially limit the conductive heat fiux by limiting 
(7 — l){mp/2kp)KQT^/'^ to a maximum value of 5x^/(8 where nip is the proton mass, k is 
the Boltzman constant, 6x is the smallest dimension of the grid cell, and 6t is the numerical 
time step calculated from the Courant condition based on all the other dynamical terms. 
Note, this limiting is a purely numerical consideration to avoid the extremely stringent 
time step that would be required by the Courant condition set by the thermal diffusion at 
large r where density decreases rapidly. This limiting is mainly altering thermal conduction 
in the region above r = 1.4i?0, to reduce it such that the time step set by the Courant 
condition for the thermal diffusion does not go below that set by the fastest dynamical 
speed. Since we are focusing on the thermal signatures that develop during the quasi- 
static phase over a time scale long compared to the dynamic time scale and inside the fiux 
rope below r = IARq where the thermal conduction is unaltered, we do not expect this 
limiting of thermal conduction to significantly alter the qualitative properties of the thermal 
signatures. Thus with the present treatment of the thermodynamics, we aim to qualitatively 
identify regions of significant heating due to the formation and dissipation of current layers 
and the resulting regions of enhanced temperature with the redistribution of heat by the 
field aligned thermal conduction. 

We solve the above compressible MHD equations using a numerical code which has its 
origin from the ZEUS3D code of Stone and Norman ( 1992a, b), but has been substantially 
modified and rewritten, in particular in the schemes used to solve the continuity, momentum 
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and the energy equations, and in the scheme for time stepping (Fan 
ease of future reference, we caU this code 



20091 [FaEl[20TT|). For 



"MFE" standing for "Magnetic Flux Eruption" , 
since it has been mainly used for modeling flux emergence from the interior into the solar 



atmosphere and corona in Cartesian geometry ( Fan |[2009 ), and for modeling CME initiation 
in the corona in spherical geometry ( Cottaar and Fan~||2009 ; Fan||2010 , 2011). Here we sum- 
marize the numerical schemes used by the MFE code for solving the above MHD equations. 
The dependent fleld variables are discretized using a staggered mesh in the spherical domain 
with r, 0, (j) coordinates ( Stone and Norman||1992a| ). A modifled, second order accurate Lax- 
Friedrichs scheme described in Rempel, Schiissler, and Knolker (2009, see eq. (A3) in that 



paper) is applied for evaluating the fluxes in the continuity and the energy equations. The 
standard second order Lax-Friedrichs scheme is used for evaluating the fluxes in the momen- 



tum equation. A method of characteristics that is upwind in the Alfven waves (Stone and 



Norman|1992b ) is used for evaluating the v x B term in the induction equation, together with 
the constrained transport scheme to ensure V-B = to the machine precision. For temporal 
discretization, the code uses a second-order accurate predictor-corrector time stepping. The 
thermal conduction term is treated via an operator splitting procedure as follows. All the 
other equations and the total energy equation (eq. [5]) excluding the thermal conduction 
term are flrst advanced one full time step (with the second-order predictor-corrector time 
stepping). Then a separate step (with a second-order predictor-corrector time stepping) is 
carried out to further update the internal energy for the thermal conduction term and adjust 
the total energy accordingly. 



The initial set up of the simulation is the same as that of Fan (2010). The spherical do 



main (see Figure[T]) representing the corona is given by r G [i?0, 5.496i?0], 9 G [57r/12, 77r/12], 
(j) G [— 7r/9.6, 7r/9.6], where Rq is the solar radius. It is resolved by a grid of 432 x 192 x 240, 
which is uniform in 6 and 0, and non-uniform in r: in the range from r = Rq to r = 1.788i?0, 
the grid size is dr = 
reaching about dr 



0.0027271 Rq = 1.898 Mm, and dr increases gradually for r > 1.788 i?©, 
= 0.09316 i?0 at the outer boundary. 



Initially, the domain is assumed to be in hydrostatic equilibrium at a uniform temper- 
ature of To = 1 MK with the density and pressure given by: 



p = po exp 
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1 



(10) 
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where po = 8.365 x 10~^^ g cm~^ is the initial density at the coronal base, and Hpo = 
{RTq/ iLi){GM(7)/ R'^)~'^ denotes the initial pressure scale height, which is about 60 Mm. The 
initial atmosphere contains a pre-existing potential arcade fleld, whose normal fleld 3^(0^9) 
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at the lower boundary (see the gray scale image on the sphere in Figure [T| is given by 
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(13) 

in which 9a = 0.05, 9t = 0.0432, 9p = 7i/2 - 9t - ^„/2, and Bq = 20G is the peak field 
strength in the arcade field. Thus the peak Alfven speed at the foot point of the arcade field 
is vao = Bo/^/4¥pQ = 1951 km s~-^ which is more than a factor of 10 greater than the initial 
sound speed of 135 km s~^. 



As was in [Fan (2010), we impose (kinematically) at the lower boundary of the domain 
(at r = Rq) the emergence of a twisted magnetic torus Btube? by specifying a time dependent 
transverse electric field E±\r=RQ that corresponds to the upward advection of the fiux tube 
at a velocity vq: 



E, 



\r=Rp 



r X 



X Btube(^©,^,0,^) I X r 



(14) 



Here the imposed velocity field on the lower boundary is a constant Vq in the area where 
the emerging tube intersects the lower boundary and zero in the rest of the area. The field 
structure Btube used for specifying E^l^^^^^ is an axisymmetric torus defined in its own local 
spherical polar coordinate system (r^, 0', (/)') whose origin is located at r = Fq = (rc^o^^o) 
of the sun's spherical coordinate system and whose polar axis is parallel to the polar axis of 
the sun's spherical coordinate system: 



Btube = V X 



r' sin 9' 



B^,{r\9') 



where 



A{T',9') = \qa^Bte^^^ t^'(r',^') 



5,'(r',^') = ^^exp (- 
r' sm 0' \ 
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(16) 
(17) 



- 7- 



where a is the minor radius of the torus, w = (r^^ + i?'^ — 2r'R' sin O^Y^'^ is the distance to the 
curved axis of the torus, with being the major radius of the torus, q denotes the angle (in 
rad) of field line rotation about the axis over a distance a along the axis, Bta/R' gives the 
field strength at the curved axis of the torus. Here we have a = O.O4314i?0, R' = O.25i?0, 
q/a = —0.0166 rad Mm~^, Bfa/R = 2.245o. The magnetic field Btube is truncated to zero 
outside of the fiux surface whose distance to the torus axis is zu = a. For specifying the fiux 
emergence via E^l^^^^^ given by (14), it is assumed that the torus' center is initially located 
at To = (ro = O.7O7i?0, = 7r/2, 00 = 0) (thus the torus is initially entirely below the 
surface), and it moves bodily towards the lower boundary at a constant velocity Vq = '^0^0, 
with Vq = 0.001'^ AO, until a time tgtp, when the emergence is stopped and E^l^^^^^ is set to 
zero. During fiux emergence, we assume that the density inside the emerging fiux tube is 
uniformly po (same as the initial density at the bottom of the domain). Thus as the tube 
is being transported into the domain via the electric field at the lower boundary, there is 
also an infiow of mass fiux po'^Or through the lower boundary in the area where the emerging 
tube intersects the boundary. When the emergence is stopped, both the electric field and the 
velocity at the lower boundary are set to zero, with no infiows or outfiows and with field line 
foot points rigidly anchored. We assume perfectly conducting walls for the side boundaries 
of the simulation domain. For the outer boundary, we use a simple outward extrapolating 
boundary condition that allows plasma and magnetic field to fiow through. These boundary 



conditions described above are the same as those used in Fan (2010) 



For the present simulation, we drive the emergence of the torus via the electric field at 



the lower boundary until tgtp = QORq/vao^ the same as case e4 in Fan (2010). However the 



difference between this simulation and those in Fan (2010) is that here we no longer assume 



the simple isothermal evolution, but instead we solve the total energy equation for a perfect 
gas with an adiabatic index 7 = 1.1 for the coronal plasma. Thus we take into account 
non-adiabatic eflFects due to heating produced by (numerical) dissipation of kinetic and mag- 
netic energies. We also incorporate thermal conduction along the magnetic field lines, as is 
described earlier in ^ Thus temperature can vary due to adiabatic expansion/compression, 
as well as due to the non-adiabatic eflFects: heating and thermal conduction. This enables 
us to identify sites of significant heating due to current sheet formation and study the resul- 
tant temperature distribution (qualitatively) during the evolution of the 3D fiux rope in the 



corona. 
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3. Results 

Figure [2] shows snapshots of the 3D evolution of the coronal magnetic field after the 
imposed fiux emergence at the lower boundary has stopped at tgtp = QORq/vao = 8.92 hour. 
Figure [3] shows the temporal evolution of the rise speed at the apex of the tracked axis of the 
emerged fiux rope (upper panel) and the total magnetic energy in the coronal domain (lower 
panel) after tgtp. The way we track the axis of the emerged fiux rope is the same as that 



described in Fan (2010, see discussion under Figure 3 in that paper). We find that the fiux 
rope first settles into a phase of quasi-static rise, and then begins to accelerate and erupt 
when a critical height of roughly 1.25/?© is reached by the tube axis. During the quasi-static 
phase after the emergence is stopped, there is no more input of Poynting fiux from the lower 
boundary and the magnetic energy in the domain is slowly decreasing (see Figure [3]) due to 



tether cutting reconnections in the current sheets that develop, as we will discussed later. 



Compared to the corresponding case e4 in Fan (2010) with the same imposed lower boundary 



conditions, the present case with a different treatment of the thermodynamics is found to 
undergo a longer quasi-static rise phase before the fiux rope reaches roughly the same critical 
height for the onset of the torus instability. As a result, at the time of the onset of eruption, 
the (free) magnetic energy is lower and the Alfven speed at the rope axis is also lower (by 
about 15%) in the present case. Thus it is found that the eruption is slightly less energetic 
with the fiux rope accelerating to a lower peak speed of about 371km/s for the eruption. 



compared to the peak speed of about 429km/s in the corresponding e4 case in Fan (2010). 
Here we focus on studying the development of the thermal signatures in the fiux rope as a 
result of current sheet formation and tether-cutting reconnections during the quasi-static rise 
phase. The tether-cutting reconnections after the emergence is stopped continue to build up 



the twisted fiux of the fiux rope (as computed in equations [14] in Fan (2010)), even no more 
Poynting fiux is coming through the lower boundary and the magnetic energy is declining. 
They enable the fiux rope to reach the critical height and eventually erupt. Therefore the 
thermal structures that develop due to the reconnections may be important signatures for 
indicating the readiness of the fiux rope to erupt, and identifying such thermal signatures in 
the observables is useful. 



Figure [4] shows snapshots of the density (in log scale) , current density normalized by 
magnetic field strength {J/B = | V x B|/5), and temperature in the central meridional cross- 
section of the fiux rope shown in the corresponding snapshots in Figure [2j In the density 
cross-section during the quasi-static rise phase (panels in the left column), the outermost low 
density layer corresponds to the pre-existing arcade fields, and the low density is produced 
by the stretching and expansion of the pre-existing arcade fields immediately surrounding 
the emerged fiux rope fields. Inside the outermost low density layer is a layer of relatively 
high density and high J/ B (indicated by the yellow arrows in Figure [i]) , which corresponds 
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to the outer boundary of the emerged flux rope. Inside the flux rope, there is a central 
low density void (indicated by the green arrows in Figure ^ which grows in size during 
the quasi-static phase. The central void is closed at the bottom with a sharp "U" shaped 
lower boundary, and is enclosed in a high density shell. The sharp lower boundary of the 
void corresponds to a "U" shaped current concentration of high J/B which extends from a 
strong central vertical current layer below (see the middle and bottom snapshots in the 2nd 
column of Figure [i]). In other words, the central void is shaped like a "lollipop" on top of 
a vertical "stick" corresponding to the central vertical current layer, and the current layer 
extends upward to surround the two sides of the void. The middle and lower images in the 
right column of Figure [4] show that the central growing void is the hottest part in the central 
meridional cross-section, hotter than the plasma in the strongest vertical current layer below 
it and the current layer surrounds it. Furthermore, this hot void is growing below the apex 
of the emerged flux rope axis marked with the yellow "x" point in the cross-sections shown 
in Figure [4j 



Figure [5] shows the 3D coronal fleld lines and the isosurface of J/5 at a value of 1/(10 dr) 
where dr is the smallest grid size in r. We flnd that at t = 9.91 hour, corresponding to the 
top row of Figure [4| the strongest portions of the current layers are two "J" shaped layers 
curving about the two legs of the flux rope at the boundaries between the flux rope legs 
and the arcade with the opposite polarity foot points. This is why the central meridional 
cross-section of J/ B (the top panel in the middle column of Figure ^ shows a double legged 
structure, with each leg belonging to one of the "J" -layers. When the "lollipop" shaped 
central hot void has developed in the central meridional cross-section (as is shown in the 
middle and bottom rows of Figure [3] for t = 18.24 hour and t = 23.79 hour), the strongest 
portion of the current layer has become a central inverse S-shaped vertical layer underlying 
the flux rope (see the middle and right panels of Figure [5| and the associated movie in the 
electronic version). The inverse-S shape is consistent with the left-handed twist of the flux 
rope (e.g. iFan and Gibson 2004). A forward-S shape for the current layer would result for 



a right-hand twisted flux rope. 

We flnd that these current layers in the simulation (as shown in Figure [5]) develop along 
topological structures identifled as Quasi Separatrix Layers (QSLs), which are regions of the 
magnetic volume where the fleld line connectivity experiences drastic variations (]Demoulin 



et al. 


1996a|bi 


Titov et al. 


2002) 



at which the fleld line linkage is discontinuous. Similar to a separatrix, a QSL divides the 
coronal domain into quasi-connectivity domains, and due to the drastic change of the fleld 
line connectivity, it is a site along which thin and intense current sheets tend to build up 
because of the distinct dynamic behaviors and evolutions of the fleld lines in each of these 
subdomains in the line-tied corona (e.g. Demoulin et al.|1996b ; Aulanier et al.|2005 ; Savcheva 
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et al. 2012b). QSLs are located by estimating the so-called squashing degree, which is 



defined as the square of the norm of the Jacobian matrix of the field line mapping from one 
foot point to the other, divided by the absolute value of the determinant of the Jacobian 



matrix (Titov et al. 2002; Titov 2007; Pariat and Demouhn 2012). QSLs correspond 



to regions of very large Q, larger than the averaged values of the rest of the domain by 



many orders of magnitude (e.g. Pariat and Demoulin 2012; Savcheva et al. 2012b). We 
have evaluated the squashing degree Q in our simulated coronal magnetic field using the 
formulation described in Pariat and Demoulin (2012, method 3 in the paper). The left 
panels in Figure [6] show the computed Q in respectively the central meridional cross-section 
(upper panel) and the cross-section at the height of r = l.O48i?0 (lower panel) in the 3D 
magnetic field shown in the right panel of Figure [5]). Compared to the corresponding 2D cuts 
of J/B shown in the right panels of Figure [6| we can see that the thin, most intense current 
layers indeed form along the QSLs with the largest Q values, which are orders of magnitude 
greater than the averaged Q in the rest of the domain. The QSLs with the highest Q values 
shown in the central meridional cross-section (top left panel of Figure [6]) correspond to the 
mid cross-section of the so-called Hyperbolic Flux Tube (HFT), a generalization of the X-line 
configuration, which divides the magnetic volume into 4 distinct domains of magnetic field 
connectivity (e.g. [Titov ||2007t [Aulanier et aL]|2005t [Savcheva etaL]|2012bD . The main central 
vertical current layer shown in the right panel of Figure [5| (whose cross-sections are shown 
in the 2D cuts in Figure |6]) , forms along the HFT and is likely a thin current sheet that can 
lead to significant reconnection even under the realistic high Lundquist number condition of 
the solar corona. 

To understand the 3-dimensional structure of the growing central hot void on top of 
the central vertical current layer seen in the cross-sections in Figure [4| we show in the top 
panels of Figure [7| the horizontal cross-sections of density and temperature at r = 1.15Rq for 
t = 23.79 hour (corresponding to the height indicated by the green arrows in the bottom row 
of Figure [4]). We see a hot channel of inverse-S shape, with the main middle segment of the 
channel tilted away from the east-west direction (or the direction of the emerging fiux rope 
axis) clock- wise by roughly 40°. There is a good anti-correlation between the temperature 
and the density of the channel, i.e. it is a hot, low-density channel. In the lower panels 
of Figure [7[ we also show the 3D morphology of the hot channel as outlined by the pink 
isosurface of temperature (at 1.2 MK) in relation to the current layers (orange isosurfaces of 
J/B) and the magnetic field lines. The hot channel is on top of the central vertical current 
layer, and the current layer extends upward to surround the two sides of the hot channel in 
the middle portion (see the movie associated with Figure [T] in the electronic version which 
show the rotating view of the 3D structures). The hot channel also threads under the axial 
field line of the fiux rope (as indicate in the bottom panels of Figure [4]). 
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If this flux rope is viewed above the hmb, with a hne of sight that is ahgned with 
the main middle segment of the hot channel to a suitable angle, then one would be able 
to observe the central void of high temperature inside the flux rope. Figure [8] shows the 
modeled images above the limb of the emission measure (panel b), line-of-sight averaged 
temperature (panel c), and SDO/AIA intensities for the 171 A (panel d), 193 A (panel e), 
and 211 A (f) fllters, by integrating through the simulation domain along the line-of-sight 
that is tilted 40° clock-wise from the east-west direction (or the direction of the axis of the 
flux rope). Also shown in the flgure is the view above the limb of the 3D fleld lines and the 
iso-surface of J/B (at a value of 1/10 rfr with dr being the smallest radial grid size) along 
the same line of sight (Figure [8^). For computing the modeled emission above the limb at 
the AIA 171 A, 193 A, and 211 A wavelength channels, we use the temperature and density 
from the simulation and carry out the following integration along the line-of-sight / through 
the domain: 

Ii = I nl{l)MT{l),n,{l))dl, (18) 



where the subscript i denotes one of the AIA wavelength channels, li denotes the intensity 
at each pixel in units of DN/s/pixel, is the electron number density, and the function 
fi{T^ne) takes into account the atomic physics and the properties of the AIA fllters. We 
have ignored the weak dependence of fi on rig, and obtained the temperature dependent 
function fi{T) using the SolarSoft routine get_aia_respoiise .pro. 

The resulting modeled emission measure and the AIA intensity images in Figure [8] show 
a central tear shaped void with a relatively sharp "U" shaped lower boundary, surrounded 
by a relatively dense shell. The central tear shaped void is the central low density, hot 
channel viewed along its main part of the length (see Figure [T]) , and the "U" shaped dense 
shell enclosing the void is produced by the upward extensions of the vertical current layer 
(see the orange isosurfaces in Figures |8^) which has formed along the HFT. Note that the 
tear-shaped void enclosed by a "U" shaped dense shell on top of a vertical dense structure is 
inside the larger flux rope structure. Such a central density sub-structure inside the flux rope 



has been observed in some of the quiescent prominence cavities at the limb (e.g. iRegnieret 



alJ2011; Reeves et al. 2012; Berger 2012). Figured shows a prominence cavity observed 



by SDO/AIA on 12 June 2012 at about 23:01 UT on the north-west limb. This is also 



the case studied by Regnier et al. (2011, see also Figure 1 in that paper for the sharper 
negative images). The AIA 171, 193, 211 images all show an ellipse shaped central void 
with a sharp "U" shaped lower boundary situated on top of the vertical prominence. The 
void is enclosed by a "U" shaped dense shell extending upward from the prominence. These 
qualitative features are similar to those seen in the modeled AIA images shown in Figure [8| 
provided that the observed prominence material is co-spatial with the vertical current layer 
in the model (Figure Isk). Prominence condensation may develop inside the denser current 



12 



layer due to radiative cooling, which is not included in our numerical simulation. The "U" 
shaped dense shell extending upward from the prominence seen most prominently in the AIA 



171 observations is also called "horns" (Berger 2012) and is found to be common features of 



quiescent prominences observed with the AIA. Our model suggests that this observed density 
feature of a dense prominence column with upward extending "horns" corresponds to the 
current sheet that forms along the HFT, when viewed from a suitable line of sight angle, as 



illustrated by Figure 10, which can be compared qualitatively with the the prominence and 



"horn" morphology shown in Figure 3 of Berger (2012). In other words, the development of 



a dense prominence column with upward extending "horns" enclosing a central void (such 
as the cases shown in [Regnier et al. (2011) and Berger| ( |2012 )) signifies the development of 
the HFT topology (and the associated tether cutting reconnections) in the cavity flux rope. 
Figure [TO] also shows two fleld lines drawn from two points in close proximity on opposite 
sides of one of the upward extensions of the central current layer. It illustrates the drastic 
difference in the connectivities of these two fleld lines which come to close proximity on 
opposite sides of the current layer extension (which is a QSL) with very different fleld line 
orientations. A 3D rotating view of the current layer and the fleld lines is available as a 



movie associated with Figure 10 in the online version 



The anti-correlation between the emission measure and the temperature for the central 
void structure (Figures [8](b)(c) ) is also consistent with the recent study of internal thermal 
properties of coronal cavities by [Reeves et al. (2012). They found a central core with lower 
density but higher temperature inside the larger coronal cavity. This central hot core is also 



consistent with the X-ray emitting core sometimes observed in coronal cavities (e.g. Hudson 



et al. 1999; Reeves et al. 2012). 



Our model suggests that this sub-structure of the hot, low density channel inside the 
larger flux rope corresponds to the reconnected, twisted flux being added to the flux rope 
via tether-cutting reconnections in the current layers which have formed along the HFT 
during the quasi-static stage (see orange surfaces in Figure [5]). Figure [TT](b) illustrates 
such a reconnected fleld line through the central hot void (see the red fleld line), which at 
t = 23.79 hour has become a fleld line that twists about the tracked axis (black fleld line) 
of the flux rope. Yet, this red fleld line has its foot points rooted in the arcade regions. 
At an earlier time t = 9.91 hour (see Figure [TT](a)), the fleld lines drawn from these foot 
points are two simple arcade fleld lines (the red fleld lines in Figure [TT](a)). Thus the 
two simple arcade fleld lines at t = 9.91 hour have undergone multiple reconnections with 
the flux rope flelds, and by time t = 23.79 hour have turn into a fleld line rooted in the 
arcade foot points but twists about the axis of the flux rope. Such addition of the twisted 
flux to the flux rope, threading under the axis, produces the expanding, low density hot 
channel in the flux rope. When viewed in a line of sight that is roughly aligned with the 
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channel, the resulting central void is found to enlarge and rise up during the quasi-static 
stage (see Figure 12 and the associated movie in the online version of the paper). As the 
critical height is reached, the central void accelerates rapidly upward and is ejected, leaving 
behind a cusp shaped brightening (see the last row in Figure 12). The development of 
the long bright current sheet trailing the ejected central void (see the movie), and the rise 
of the cusp shaped brightening are more prominent in the 193 and 211 channels due to 
their sensitivity to higher temperatures. This evolution sequence can be compared with the 
observed evolution of the prominence cavity studied in Regnier et al. (2011), which begins 
to erupt at about 03:24 UT on June 13 2010. A composite movie of the AIA observation 
in the 171, 193 and 211 channels of the evolution of this cavity on the north-west limb 
during the period from 18:00 UT 12 June 2010 to 21:30 UT 13 June 2010 can be viewed 
at http: //sdowww. Imsal . com/suiitoday/index.html?suntoday_date=2010-06-12#. The 
observed evolution show many qualitative similarities with the modeled one, including the 
morphology of the central void enclosed in the U-shaped dense shell on top of a dense column 
during the quasi-static phase, and the ejection of the void and the post-eruption brightening. 



4. Conclusions 



We have carried out an MHD simulation of the quasi-static rise and the onset of eruption 
of an anchored coronal flux rope. Earlier in the simulation, the emergence of a twisted flux 
rope is imposed at the lower boundary into a pre-existing coronal potential arcade fleld. 
Then the emergence is stopped at the lower boundary and the coronal magnetic fleld is 
allowed to settle into a quasi-static equilibrium. It is found that the flux rope flrst settles 
into a quasi-static rise phase with an underlying sigmoid-shaped current layer developing 
along the HFT, which is argued to be a likely site for the formation of thin current sheets 



and magnetic reconnections (e.g. Titov 2007; Aulanier et al. 2005; Savcheva et al. 2012b). 



Tether cutting reconnections in the current layer effectively add twisted flux to the flux rope 
and allow it to rise quasi-statically (even as the free magnetic energy is decreasing) and the 
flux rope undergoes dynamic eruption when it reaches the critical height for the onset of the 



torus instability (Aulanier et al. 2010; Fan 2010). In the present simulation we study the 



thermal signatures that result from the development of the HFT topology and the associated 
tether cutting reconnections in the the current sheets that form along the HFT. We solve 
the total energy equation in conservative form for a perfect gas with an adiabatic index 
7 = 1.1, and thus incorporate the non- adiabatic heating resulting from the dissipation of 
kinetic and magnetic energies due to the formation of intense current layers. Furthermore, 
thermal conduction along magnetic fleld lines is also included in the energy equation so as 
to incoporate the redistribution of heat along the fleld lines. In this way, we are able to 
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identify qualitatively the thermal signatures in the flux rope due to heating produced by 
the tether-cutting reconnections in the current layers that form during the quasi-static rise 
phase of the flux rope 

Intense current layers are found to develop during the quasi-static phase along QSLs, 
which are topological boundaries where the magnetic fleld connectivity undergoes drastic 
variations and are likely sites for the formation of thin current sheets and signiflcant magnetic 
reconnections (Demoulin et al. 1996a||b; Titov et al. 2002). This has been shown also in 



previous MHD simulations of coronal flux rope evolution as well as nonlinear force-free fleld 
models (NLFFF) of observed sigmoid active regions (e.g. Aulanier et al.|2005 , 2010; Savcheva 
et al.| 20^12a|[b ). The strongest portions of the current layers that form are initially two "J" 
shaped layers curving about the two legs of the flux rope at the boundaries between the 
flux rope legs and the arcade with the opposite polarity foot points (see the left panel of 
Figure [5]). Later the strongest portion of the current layer becomes a central inverse S- 
shaped vertical layer that forms along the HFT (see right panel of Figure [5] and the 2D 
cuts shown in Figure |6]) . Reconnections in the current layers are found to add to the flux 
rope twisted flux threading under the axis, creating an expanding, low density hot channel 
in the flux rope above the central vertical current layer. When viewed along a line of sight 
that is roughly aligned with the channel, it manifests as a central hot void on top of a 
relatively dense structure corresponding to the central warped, vertical current layer, which 
also extends upward on two sides of the void (see the J /B image at t = 18.24 hour in Figure 
4]), producing the "U" shaped dense shell (or "horns") enclosing the void in the modeled EUV 
images (see e.g. the t = 12.89 hour and t = 18.34 hour panels of the modeled AIA 171 images 
in Figures [T2| ) . The dense vertical column and the upward extended horns are essentially the 
manifestation of the HFT (see the 2D cut in the top left panel of Figure |6]) . Such thermal 
structures inside the flux rope are consistent with several AIA EUV observations of coronal 
prominence cavities where it is found that an elevated, ellipse or tear-shaped central void is 
situated on top of the vertical column of dense prominence, and the void is enclosed by a 
"U" shaped dense shell (sometimes called "horns") extending upward from the prominence 



column (e.g. [Regnier et al.|2011; Berger 2012). The central hot channel produced by the 



tether cutting reconnections can also explain the observational result of a central hot core 



with low emission measure inside a coronal cavity by Reeves et al. (2012). 



Although our numerical simulation does not include radiative cooling and thus cannot 
model the process of prominence formation, it is likely that the prominence plasma is co- 



spatial with the dense current layer (e.g. van Ballegooijen and Cranmer 2010; Low et al. 



2012j). Our MHD simulation also does not model explicitly the background coronal heating 
but instead simply assume a high temperature coronal plasma with an adiabatic index of 
7=1.1, without explicitly incorporating the coronal heating and radiative cooling processes. 
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However our model does take into account localized heating due to the formation and dis- 
sipation of current sheets in the corona and the redistribution of the heat by field aligned 
conduction in the corona. The lower boundary of our simulation domain is placed at the 
coronal base instead of the chromosphere, and after the emergence is stopped, the velocity 
field is set to zero on the lower boundary and no more mass fiux is allowed through it into 
the coronal domain. This precludes the important effect of chromospheric evaporation due 
to conductive heat fiux coming down from the corona, which may significantly affect the 
density of the heated coronal fiux tubes. Thus there are great uncertainties in the density 
and temperature obtained due to our still highly simplified treatment of the thermodynam- 
ics, especially for the heated, twisted field lines in the hot channel that form as a result of the 
tether cutting reconnect ions. One main question would be whether the heated field in the 
hot channel would remain a void (i.e. of lower density) in contrast to the surrounding as is 
in our model when the effect of chromospheric evaporation produced by the downward heat 
conduction is taken into account. In regard to this question, one consideration is that the 
reconnected field lines in the hot channel are generally long field lines (see e.g. the red field 
lines in panel [lT](b) and Figure 10). The scaling law derived based on the static, thermal 



equilibrium coronal loop model of Rosner et al. ( 1978| ) gives Tlnax 1.4 X 10^ (pL)^/^ where 



^max, P and L denote respectively the maximum loop temperature, the loop pressure, and 
the loop length. This would mean that the plasma density of the loop p oc T^^-^L~^ . Thus 
even though the field lines in the hot core are heated to a mildly increased T^ax^ to roughly 
1.2 MK from the 1 MK original coronal temperature, these field lines are on the other hand 
considerably longer in L compared to the field lines immediately outside the hot core: for 
example, in Figure [TOj the red field line inside the central void is about 1.9 times the blue 
field line immediately outside of the horn. Therefore the above scaling law of static thermal 
equilibrium coronal loops would indicate that the effect of the longer length outweighs the ef- 
fect of the high temperature, and the resulting density inside the central hot channel remains 
smaller compared to the surrounding, even with the effect of the chromospheric evaporation 
included. 

Our 3D MHD simulations with the simplified thermodynamics is an initial step to under- 
stand qualitiatively the immediate thermal effect of the formation and dissipation of current 
sheets along the HFT in the corona. Our results suggest that the observed feature of promi- 
nence with upward extending "horns" enclosing an elevated central void seen in AIA 171 
observations of coronal cavities is a signature of the development of the HFT topology and 
the associated tether-cutting reconnections in the current sheets that form along the HFT. 
It is likely that the current sheet forming along the HFT remains the region of the highest 
density because of the significant local compression by the Lorentz force even when more 
realistic treatment of the thermodynamics are included, and our results about the thermal 
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sigature, i.e. relatively hot, low density central void on top of and surrounded by the higher 
density current sheet formed along the HFT, remain qualitatively true. Of course 3D MHD 
simulations with a full treatment of the thermal dynamics that incorporates explicit coronal 
heating, radiative cooling, as well as field aligned thermal conduction, and with the lower 
boundary set at the chromosphere, are needed to confirm the findings here and provide quan- 
titative determination of the density and temperature of the coronal plasma. As has been 
found in ID hydrodynamic simulations with a full treatment of the thermodynamics along 



long dipped coronal loops (e.g. Karpen and Antiochos 2008), complex, thermal nonequi- 



librium formation and dynamic evolution of prominence condensations can develop along 
the coronal loops depending on the spatial and temporal properties of the heating profiles 
prescribed. 

Here we have presented the MHD simulation of the quasi-static evolution and onset of 
eruption of a particular coronal fiux rope configuration. However the results presented here 
of the development of the HFT topology, and the associated current sheet formation and 
tether cutting reconnections along the HFT are quite general results for the evolution of 
line-tied coronal fiux ropes as has been demonstrated by several previous MHD simulations 



as well as NLFFF models of sigmoid active regions based on observations (e.g. Aulanier et 



al. 2005, 2010; Savcheva et al. 2012a|[b ). These models suggest that the development and 



the elevation of the HFT may be an indication of the extent to which the tether-cutting 
reconnections have progressed and hence the readiness of the fiux rope to erupt. 
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Fig. 1. — The initial configuration of the simulation, same as that used in Fan (2010), see 
text for details. 
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Fig. 2. — Snapshots of the 3D coronal magnetic field lines after the imposed fiux emergence 
is stopped. A movie of the magnetic field evolution together with the evolution shown in 
Figure [4] below is also available in the electronic version of the paper 
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Fig. 3. — Rise velocity at the apex of the tracked axis of the emerged flux rope (upper panel) 
and the total magnetic energy in the simulation domain (lower panel) as a function of 
time after the imposed flux emergence has stopped. The dotted line marks the magnetic 
energy of the corresponding potential field, Ep^ which remains unchaged after the emergence 
has stopped. The eruption causes a roughly 43% reduction of the free magnetic energy. 
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Fig. 4. — Snapshots of the density (in log scale), current density normalized by magnetic 
field strength {J/ B = |V x B|/5), and temperature in the central meridional cross-section 
of the flux rope shown in the corresponding snapshots in Figure [2j Here the density is in 
units of po = 8.365 x 10"-^^, and J/B is in units of I/Rq. The features indicated by the 
green and the yellow arrows and the yellow 'X' points are discussed in the text. A movie 
of the evolution shown in this flgure together with the evolution shown in Figure [2] is also 
available in the electronic version of the paper 
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Fig. 5. — The 3D coronal field lines and the iso-surface of J /B dX value of 1/(10 dr) where 
dr is the smallest grid size in r at times t = 9.91 hour, t = 18.24 hour, and t = 23.79 
hour. An animation of the 3D structures with a rotating perspective is also available in the 
electronic version of the paper. 
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Fig. 6. — Estimation of the squashing degree Q in respectively the central meridional cross- 
section (upper left panel) and the cross-section at the height of r = l.O48i?0 (lower panel) in 
the 3D magnetic field shown in the right panel of Figure [5} and the corresponding 2D cuts 
of J/B (right panels). 
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Fig. 7. — Top panels show the horizontal cross-sections of density and temperature at r = 
1.15i?0 for t = 23.79 hour, corresponding to the height indicated by the green arrows in the 
bottom row of Figure [4j The lower panels show two perspective views of the 3D magnetic 
field lines, the pink iso-surface of temperuature at 1.2 MK, which outlines the central hot 
channel, and the current layers outlined by the orange iso-surfaces of J/B at the level of 
l/(10rfr), with dr being the smallest radial grid size. A movie showing the rotating view of 
the 3D structures in the lower panels is available in the electronic version of the paper. 
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Fig. 8. — (a) The 3D field lines and the iso-surface of J/S at a value of 1/(10 dr) with 
dr being the smallest radial grid size, viewed above the limb along the line of sight that is 
tilted by 40° clock- wise from the east-west direction (or the direction of the axis of the fiux 
rope). Also shown are modeled images of the emission measure (b), line-of-sight averaged 
temperature (c), and SDO/AIA intensities at 171 A (d), 193 A (e), and 211 A (f) channels, 
by integrating through the simulation domain along the same line-of-sight. 




Fig. 9. — A prominence cavity observed by SDO/AIA on 12 June 2012 at about 23:01 UT 
on the north-west hmb. See also Figure 1 of Regnier et al. (2011) which shows sharper 
inverse-color images of the same cavity at a later time on 13 June 2012 at about 03:36 UT. 



- 28 - 



t = 23.79 hour 




Fig. 10. — A 3D view of the current layers as outlined by the iso-surface of J/5 at a value 
of l/(10rfr), with dr being the smallest radial grid size, and two field lines traced from two 
points in close proximity on opposite sides of one of the upward extensions of the central 
vertical current layer. A movie showing the 3D rotating view is availabe in the online version. 
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(a) t = 9.91 hour (b) t = 23.79 hour 




Fig. 11. — 3D views of selected field lines from two perspectives. Black field lines are the 
tracked axis of the fiux rope. The red field lines are field lines traced from two fixed foot 
points in the arcade region. At t = 9.91 hour, the red field lines are two simple arcade loops. 
At t = 23.79 the red field line traced from the same two foot points in the arcade region has 
become a field line that threads under and twists about the axis of the fiux rope, and it goes 
goes through the central hot void. 
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Fig. 12. — Snapshots of the modeled AIA intensity images at 171 A, 193 A, and 211 A 
channels at 3 time instances, showing the expansion and rise of the central void during the 
quasi-static phase, and its final ejection and the post eruption brightening. A movie of this 
evolution is available in the online version of the paper. 



